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Abstract: Wavelength-dependent point spread functions (PSFs) violate an implicit assumption 
in current galaxy shape measurement algorithms that deconvolve the PSF measured from stars 
(which have stellar spectral energy distributions (SEDs)) from images of galaxies (which have 
galactic SEDs). Since the absorption length of silicon depends on wavelength, CCDs are a potential 
source of PSE chromaticity. Here we develop two toy models to estimate the sensitivity of the 
cosmic shear survey from the Earge Synoptic Survey Telescope to chromatic effects in CCDs. We 
then compare these toy models to simulated estimates of PSE chromaticity derived from the ESST 
photon simulator PhoSim. We find that even though sensor contributions to PSP chromaticity are 
subdominant to atmospheric contributions, they can still significantly bias cosmic shear results if 
left uncorrected, particularly in the redder filter bands and for objects that are off-axis in the field 
of view. 
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1. Introduction 

All galaxy shape measurement algorithms for ground-based telescopes depend, in one way or an¬ 
other, on measuring the point spread function (PSF) from the images of stars to then deconvolve 
from the images of galaxies. The implicit assumption in this approach is that the PSF affecting stars 
and galaxies is locally the same. If the PSF depends on wavelength, however, then this assumption 
is violated since stars and galaxies have different spectral energy distributions (SEDs). 

Several papers have addressed aspects of this effect to date. Cypriano et al. (2010) OH] analyze 
the impact of the wavelength-dependent diffraction limit of the Euclid PSF, Plazas & Bernstein 
(2012) [0] investigate differential chromatic refraction for ground-based telescope PSFs such as 
DES and ESST, and Meyers & Burchat (2015) [0] study both differential chromatic refraction and 
chromatic seeing for ground-based telescopes. 

Chromatic effects are also likely to arise in CCDs due to the wavelength dependence of the 
absorption length in silicon [0. ESST will primarily measure galaxy shapes in the r- and /-band 
filters, for which the typical silicon absorption lengths are ~ 2-5 microns and ~ 5-20 microns, 
respectively, which can be compared to the ESST CCD thickness of 100 microns and pixel width 
of 10 microns. In this article, we investigate some consequences of the wavelength-dependent 
absorption depth on the PSE. 

In Section 0 we briefly review some analytic formulae for estimating shape measurement bi¬ 
ases due to PSE misestimates, including those arising from chromatic PSEs. We then use these an¬ 
alytic expressions to investigate two toy models of chromatic CCD PSEs in Section 0, and roughly 
estimate the tolerance of ESST weak lensing science to generic chromatic CCD effects. In Section 
0 we use the ESST photon simulator package PhoSim to estimate the magnitude of chromatic PSP 
effects for ESST and present our conclusions in Section 0. 
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2. Analytics 


Weak gravitational leasing is frequently analyzed through its effect on combinations of the second 
central moments of a galaxy’s surface brightness distribution. For a detailed description of this 
process, see Ref. [01. We briefly summarize the salient points for the present study here. 

The zeroth, first, and second unweighted central moments of the surface brightness profile 
I{x,y) are 

f = I I{x,y)dxdy, (2.1) 

l^ = jJ I{x,y)lJ-dxdy, (2.2) 

and 

I^v = yjl{x,y){l^-fl){v-v)dxdy, (2.3) 

where jj. and V refer to either x or y, p. and V indicate either the x or y centroids of the profile, 

and / is the total flux of the profile. In practice, when applying 0-0 to real imaging data, a 
weight function w(x,y) is either explicitly (for moment-measuring algorithms) or implicitly (for 
model-fitting algorithms) included as a factor in each of the integrands to prevent noise in pixels 
at large radii from causing the integrals to diverge. The second central moments of interest for 
weak leasing are those of galaxies before convolution with the PSF of the atmosphere, optics, and 
sensors. The observed moments (after PSF convolution) are related to the intrinsic moments by 
which holds exactly only for unweighted second moments, but is still a useful 
approximation (correct to within a factor of a few) when using practical weight functions or model¬ 
fitting approaches. 

Important combinations of second moments for weak leasing include the second-moment 
squared radius r^, and the complex ellipticity X = ^X 2 - 


Xi 


^xx + lyy, 

^xx ^yy 
^xx + Iyy'’ 


(2.4) 

(2.5) 


Xi 


21 

jLLxy 

^XX + lyy 


(2.6) 


Under the assumption that galaxy ellipticities are intrinsically isotropically distributed, the cos¬ 
mological signal of interest, the reduced shear g, is related to the ensemble average ellipticity by 

{x^‘^^)-2g. 

Chromatic PSF effects are a subset of a larger class of effects in which one misestimates the 
size or shape of the PSF to be deconvolved from galaxy images. Mathematically, we can describe 
this PSF misestimate via the error in the second moments: Alj^y^ = where 

are the true PSF moments appropriate for deconvolving a particular galaxy and ’ are the PSF 
moments estimated from one or more stars. 


To characterize the impact of PSF misestimation, we follow the literature and parameterize the 
bias induced in the reduced shear via a multiplicative and additive term: 


gi=gi(l+mi)+Ci, i = l,2, 


(2.7) 
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where g is the estimator for the true reduced shear g, and we have assumed that gi {§2) is indepen¬ 
dent of g2 (gi). The results of propagating PSE misestimates into the multiplicative and additive 

shear biases are [0 


-(A/PSP + A/^f) 

mi = m 2 = 2 > 

(2.8) 

''gal 


A/PSP-A/;f 
= 2F 

^'^gal 

(2.9) 

A/PSP 


C2- f . 

(2.10) 


gal 


The typical source galaxy size for LSST is about ^gai ~ (0-3 arcsec)^ [0. 

Requirements on how well one needs to know m and c have been derived by several authors 
for surveys similar to and including LSST [0,0,0, ED- Interpreting these requirements can be tricky, 
however, and we refer the reader to Ref. [0] for more detailed discussion. Following Ref. [0, we set 
('M)req = 0.003 (that is, the requirement on our knowledge of the mean of m in each tomographic 
redshift bin), as the point at which the systematic uncertainties of LSST weak lensing become 
equal to the statistical uncertainties. Similarly, we set Var(c)jgq = 1.7 x 10^^ as a conservative 
upper bound on the point at which uncorrected additive systematic uncertainties equal statistical 
uncertainties for LSST weak lensing. 


3. Toy models 


In this section we construct a toy model for chromatic CCD PSF effects. To isolate chromatic 
CCD effects from other types of chromatic effects, in our model we will assert that the PSF con¬ 
tribution coming from the atmosphere and telescope optics is achromatic. For LSST, the atmo- 
spheric-i-optical PSF is dominated by the atmosphere, which we will model as a Moffat profile 
with FWHM of 0.65-arcsec and j8 parameter of 5, which corresponds to Ixx = lyy = 0.12 arcsec^. 
We model the typical contribution of the CCD to the PSF as a Gaussian with FWHM of 0.19 arc- 
sec, which is appropriate for charge diffusion in the LSST CCDs 0, 0,110, and corresponds to 
Ixx = lyy = 0.0065 arcsec^. 

Mathematically, our model is 


tPSF / t \ TPSF,telescope+optics 


,PSF,CCD 


(A). 


(3.1) 


For a given SED, we can compute an effective PSF as 

/PMoc |/;;f(A)R(A)S(A)AdA, (3.2) 

where R(A) is the total system throughput, S(A) is the SED, and the normalization is set such that 
the integral of the effective PSE over all space is 1. Eq. can then be used to compute the PSE 
misestimate and resulting multiplicative or additive shear bias for a given pair of stellar and galactic 
SEDs. 


- 3 - 





3.1 Wavelength-dependent CCD PSF size 


To make the CCD PSF ehromatie, we will slowly vary the seeond moments about their fidueial 
values as a funetion of wavelength. Our first toy model varies the size of the CCD PSF as a 
funetion of wavelength: 

(A) = /Jf(Ao) + ^+(A - Ao), (3.3) 

/Jf(A) = ’''^(Ao) + ^+(A - Ao). (3.4) 

We ehoose Ao = 700 nm to fall in between the LSST r- and /-bands. To quantify the sensi¬ 
tivity of LSST shape measurement to a wavelength-dependent CCD PSF size, we ehoose j8+ = 
10^^ aresee^/nm, whieh produees a CCD PSF with FWHM about 10% smaller at the blue edge of 
the /-band filter than at the red edge. This roughly eorresponds to d(FWHMccD)/ dA ~ 0.14 mas/nm, 
whieh, as we will see in Seetion 4, is about the right size for the LSST /-band. 

Sinee this model ehanges only the sum Ixx + lyy, but leaves the differenee Ixx — lyy and Ixy 
fixed, if generafes only a mulfiplieafive bias. Figure 1 shows fhe magnifude of Ibis mulfiplieafive 
bias relafive fo a G5v sfar for a variefy of sfellar and galaefie SEDs. The shaded region in fhe 
plof indieafes fhe LSST requiremenf for how well we need fo know fhe average mulfiplieafive 
bias in eaeh fomographie redshiff bin so fhaf sysfemafie uneerfainfies remain af or below sfafisfieal 
uneerfainfies. Sinee fhe differenl SEDs roughly span Ibis requiremenf range, we eonelude fhaf 
j3+ ~ 10^^ aresee^/nm or d(EWHMccD)/dA ~ 0.14 mas/nm deseribes fhe poinf af whieh fhe 
uneerfainfy from Ibis sysfemafie effeef rivals fhe sfafisfieal uneerfainfy of ESST weak lensing. 

3.2 Wavelength-dependent CCD PSF ellipticity 

We also eonsfruef a somewhat orthogonal toy model in whieh the elliptieity, but not the size, of the 
CCD PSF is ehromatie: 


=/Jf’CCD(A o) + ^_(^ _Ao), (3.5) 

/Jf’^®(A) ’CC°(Ao) -^-(A - Ao). (3.6) 

Using the same Ao = 700 nm as before and setting j8_ = 10^^ aresee^/nm results in a CCD PSF 
that ehanges in elliptieity from about —0.1 at the blue edge of the /-band to about 4-0.1 at the red 
edge of the /-band. This is roughly equivalent to d;^i ccd/ dA ~ 0.0014 nm^^, whieh, as we will 
see in Seetion 4, is a faetor of a few smaller than observed in ESST /-band simulations. 

This model induees an additive bias in the elliptieity eomponent, but not in the Xi elliptieity 
eomponent. It does not induee a multiplieative bias in either elliptieity eomponent. Figure 2 shows 
the magnitude of this additive bias relative to a G5v star for a variety of stellar and galaetie SEDs. 
The shaded region in the plot indieates the ESST suffieient-but-not-neeessary requirement for how 
well we need to eonstrain the standard deviation of the additive bias in eaeh fomographie redshift 
bin so that systematie uneerfainfies remain at or below statistieal uneertainties. In this ease, we see 
that the requirement bar is a faetor of a few smaller than the range spanned by the various SEDs. We 
therefore eonelude that j3_ ~ 3 x 10^^ aresee^/nm or equivalently d;^i ccd /dA ~ 5 x 10^^ nm^^, 
deseribes an upper limit to the point at whieh the uneerfainfy from this systematie effeet eould rival 
the statistieal uneertainty of ESST weak lensing. 


- 4 - 




redshift 


Figure 1. Multiplicative shear bias arising from a toy model with linear wavelength dependence 

= )3+, calculated for LSST i-band. Star symbols at redshift 0 represent 
stellar SEDs from [Oil . Curves represent galactic SEDs from [IT^ and lHHl . The gray band indicates how 
well one needs to know the multiplicative bias in each LSST tomographic redshift bin so that systematic 
uncertainties are smaller than statistical uncertainties. The biases and requirement are similar in the r band, 
which is the other band planned for shape measurement with LSST. 


4. PhoSim estimates of LSST CCD chromaticity 

To connect our toy models to the actual LSST survey, we use the LSST photon simulator (PhoSim) 
El, which comprises one part of the LSST simulations package. PhoSim works by drawing pho¬ 
tons from sources in a provided catalog, and then manipulating them one-by-one as they traverse 
the atmosphere, telescope and camera optics, and convert into electron-hole pairs inside the LSST 
sensors with a wavelength-dependent absorption length [1131] . These photo-electrons are then fur¬ 
ther followed as they diffuse towards the charge collection areas of each simulated pixel and are 
eventually read out by the simulated sensor electronics. In addition to its use as a high-fidelity sim¬ 
ulator of LSST before the actual telescope is available, PhoSim is a valuable resource for probing 
individual physical effects, which can be switched on and off on demand in the software. 

Here we use PhoSim to probe chromatic PSF effects in the LSST sensors. We create a series 
of catalogs, each of which consists of a 8 x 8 grid of stars with spacing 0.01 degrees. All of the 
stars in a given catalog have the same monochromatic SED - i.e., the flux is non-zero at only a 
single wavelength. This wavelength is then scanned from 325 nm to 1050 nm in steps of 25 nm 
to produce different catalogs. We tuned the amplitude of each monochromatic SED such that the 
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Figure 2. Additive shear bias arising from toy model with linear wavelength dependence d/o: ’ / dA = 

— d/yj, ’ /dA = )3_, calculated for LSST /-band. Star symbols at redshift 0 represent stellar SEDs from 
[ITTll . Curves represent galactic SEDs from [O] and flTHl . The gray band indicates the sufficient requirement 
for LSST - i.e. the point at which systematic uncertainties will definitely not exceed statistical uncertain¬ 
ties. The biases and requirement are similar in the r band, which is the only other band planned for shape 
measurement with LSST. 


resulting stellar images each contained about 20,000 counts. We replicate this pattern of stars at 
three different field angles: on axis, 0.725 degrees off axis (about half way from the center to the 
edge of the focal plane), and 1.45 degrees off axis (near the edge of the focal plane). The off- 
axis positions are azimuthally positioned along the positive “x”-direction. The zenith angle of the 
telescope boresite is always set to an altitude of 89 degrees in the simulation, which is close enough 
to the zenith that differential chr'omatic refraction is negligible across the entire field of view. 

To isolate effecls arising in fhe sensors, we run PhoSim for each inpuf cafalog in fwo differ- 
enf modes: once wifh fhe afmosphere, optics, and sensors all fumed on, and once wifh only fhe 
afmosphere and optics fumed on. To simplify our analysis, we also furn off a number of addifional 
physics effecls, including fhe sky background, safurafion, blooming, fracking errors, oplics per- 
lurbalions, clouds, quanfum-efficiency varialions, optical surface imperfeclions (dusf), diffraclion, 
large angle scattering, laferal elecfric fields (e.g., free rings and edge effecls), charge sharing (i.e., 
fhe brighter-fatter effecl), imperfecl pixel boundaries, and shutter errors. 

In fhe resulfing images, we fif a Gaussian profile wifh six free paramelers (xq, yo^ /^ Zi> X 2 , 
FWHM) lo each sfar. (We also fried tiffing a Moffal profile fo fhe sfars; fhe resulls are essenlially fhe 
same.) For each wavelengfh and mode, we store fhe mean and slandard devialion of fhe paramelers 
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Figure 3. The FWHM of the LSST CCD PSF estimated using PhoSim as a function of wavelength and for 
three different field angles. Different colors indicate different simulated LSST filters (ugrizy). 


of the 64 stars in the 8x8 grid. In Figure 3, we plot the quadrature difference of the mean FWHM 
with the sensors turned on and with the sensors turned off as a function of wavelength and for the 
three different field angles. By plotting the quadrature difference, we are isolating the contribution 
to the PSF from the CCDs at each wavelength. We see that CCD PSF FWHM is constant at bluer 
wavelengths before rising around 800 nm, roughly independent of field angle. In Table |I], we lisf 
fhe measured values of d(FWHMccD)/dA for each filler and field angle. We see lhal in /-band, Ihe 
measured slope is only aboul a factor of Iwo lo Ihree smaller lhan fhe requiremenl we eslablished 
in Section ^TTi The slope is smaller in r-band, Ihough fhe uncertainties on fhe measuremenls are 
significanl here. 

To invesligale fhe CCD PSF elliplicily, we converf our measuremenls of X\, X 2 , and FWHM 
info second cenlral momenls sublracl fhe second cenlral momenfs wilh and wilhoul fhe sen¬ 
sor fumed on, and converf fhe residuals back info Xi and Xi- Figure 4 shows fhe X\ elliplicily 
componenf againsl wavelenglh for several field angles (i.e., for several angular disfances from fhe 
center of fhe focal plane), and for an azimufhal direclion of displacemenf along fhe “h-x” direction 
of fhe focal plane. We see lhal fhe CCD PSF elliplicily is nearly consfanf on-axis, bul becomes 
progressively larger wilh wavelenglh off-axis. Al Ibis azimufhal angle, fhe Xi residual componenf 
(no! plofted) is uniformly consistenf wilh zero for any field angle. Repealing fhe above experimenl 
for azimufhal angles along fhe “x=y” direclion (“H-y” direction) we find lhal fhe Xi (A 2 ) elliplicily 
is zero and Xi (At) elliplicily is positive (negative) off axis, which slrongly implies lhal Ihe effecl 
is radially symmelric in fhe focal plane. In Table [J we lisf Ihe measured values of d;^i ccd/ dA for 
each filler and field angle for fhe case when Ihe azimufhal displacemenf is along fhe “h-x” direction. 
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Figure 4. Ellipticity Xi of the LSST CCD PSF estimated using PhoSim as a function of wavelength and for 
three different field angles. Different colors indicate different simulated LSST filters (ugrizy). 


Table 1. Estimated LSST CCD PSF chromaticity 


Filter d(FWHMccD)/d‘^ (mas/nm) d;^i.ccD/dA (/im 



o 

O 

O 

0.725° 

1.45° 

o 

O 

o 

0.725° 

1.45° 

u 

0.08 ± 0.03 

0.03 ± 0.03 

0.02 ± 0.03 

0.4 ± 0.3 

0.6 ± 0.3 

0.4 ± 0.4 

g 

0.02 ± 0.02 

0.00 ± 0.04 

0.01 ± 0.02 

0.1 ±0.2 

0.1 ±0.3 

0.2 ± 0.2 

r 

0.03 ± 0.03 

0.00 ± 0.02 

-0.02 ± 0.04 

-0.2 ±0.1 

0.2 ± 0.3 

0.2 ± 0.2 

i 

0.04 ± 0.04 

0.09 ± 0.05 

0.05 ± 0.02 

0.1 ±0.3 

0.6 ± 0.3 

1.0 ±0.4 

z 

0.65 ± 0.20 

0.71 ±0.18 

0.51 ±0.11 

0.3 ± 0.2 

0.7 ± 0.2 

1.5 ±0.2 

y 

0.94 ± 0.24 

0.96 ±0.19 

0.79 ±0.14 

-0.1 ±0.2 

-0.1 ±0.2 

0.4 ± 0.3 


While the slope on-axis is consistent with zero, off axis the slope is a factor of one to three larger 
than our requirement of ~ 0.5/im^' from Section 

5. Discussion 

We expect the response of CCDs to be chromatic due to the wavelength-dependent absorption 
length of silicon. Redder photons penetrate farther into the silicon on average before converting 
into electron-hole pairs. Naively, one might then conclude that photo-electrons generated by redder 
photons would on average have less opportunity to laterally diffuse and would ultimately produce 
a smaller PSF. However, Figure 3 suggests that the opposite is true. The naive argument above 
neglects the fact that the LSST beam is extremely fast - f/1.2. Photons arrive at the CCD with 
a mean angle of incidence around 20 degrees, with a range of R3 ±5 degrees due to the annular 
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3.2. 



















Figure 5. Cartoon figure of LSST beam and CCD geometry. Fifteen 10-micron wide and 100-micron deep 
pixels are shown in each cartoon. The f/1.2 beam is shown approaching the CCD from above for bluer filters 
(left) and redder filters (right). Blue photons convert into electron-hole pairs almost immediately, and the 
photo-electrons diffuse the entire depth of the CCD (shown in yellow). Red photons may penetrate deeper 
into the CCD bulk before converting. Combined with the diverging beam, which is still f/4.5 inside the 
silicon, the longer absorption length of redder photons leads to a larger PSF. Note that the position of best 
focus is a competition between beam divergence and charge diffusion. For bluer filters, the best focus is at 
the CCD surface, while for redder filters it is inside the bulk of the CCD 0]. 


primary mirror, which leads to a mean angle of refraction around 5 degrees inside the silicon. 
As shown in Figure 5, blue photons convert into electron-hole pairs almost immediately, and the 
photo-electrons diffuse across the entire depth of the CCD. Redder photons convert further inside 
the CCD on average, which, combined with the lateral component of their momentum, yields a 
larger PSF. Note that the position of best focus is inside the bulk of the CCD for redder filters as a 
result 0]. 

Similarly, even the chief ray of the LSST beam (though blocked by the secondary mirror and 
camera) is inclined to the CCD normal at an off-axis field position. In other words, at an off-axis 
field position, the unblocked rays do not average together to normal incidence. This is likely the 
source of the wavelength-dependent ellipticity at off-axis field positions in Figure 4. 

Although the CCD contribution to the LSST PSF is subdominant to the atmospheric contri¬ 
bution by a factor of a few in FWHM, both effects we have investigated here appear to be large 
enough to introduce significant biases in cosmic shear results if left uncorrected. Higher-order 
CCD effects, such as tree-rings, charge sharing, or edge effects may also be chromatic and must be 
investigated, especially in combination with the effect investigated here and with atmospheric and 
optical chromatic effects, since there may be non-trivial interactions among the different effects. 

If the SED of each object and the wavelength dependence of the PSF are known, then one can 
compute potential chromatic PSF misestimations and remove them from the analysis, galaxy-by¬ 
galaxy and star-by-star. Fortunately, the six-band photometry that LSST will obtain of each source 
will significantly constrain the SED across the r- and /-band filters where shape measurements are 
most effective. In Ref. |0], we have shown that a significant reduction in bias can be achieved 
by training a machine-learning algorithm to predict chromatic biases from photometry. Such a 
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correction requires an accurate and precise knowledge of the wavelength dependence of the PSF, 

however, which motivates continued study of chromatic effects both in simulation and in data where 

possible. 
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